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Abstract 

The behaviour of magnetohydrodynamic (MHD) natural 
convection flow from a porous vertical plate is studied with 
respect to variable viscosity. To model the flow, the 
governing boundary layer equations are first transformed 
into a non dimensional form. The resulting non linear system 
of partial differential equations is then solved numerically 
using finite difference method together with Keller-Box 
scheme. The study shows that the variation of viscosity 
influences the surface shear stress (skin friction coefficient) 
and the rate of heat transfer (local Nusselt number). It has an 
effect on velocity as well as temperature profiles. Studied has 
been performed with a selection of parameters set consisting 
of magnetohydrodynamic parameter M, viscosity variation 
parameter y, Prandtl number Pr. Results are shown 
graphically and tabular form and the comparison agrees well 
with a published paper. 
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Nomenclature 

a r Rossel and mean absorption co-efficient 

Cf Local skin friction coefficient 

Cp Specific heat at constant pressure 

f Dimensionless stream function 

g Acceleration due to gravity 

k Thermal conductivity 

Nux Local Nusselt number 

Pr Prandtl number 

Q Heat generation parameter 

qw Heat flux at the surface 

t/ c Conduction heat flux 

<7, Radiation heat flux 

Rrf Radiation parameter 

T Temperature of the fluid in the boundary layer 


T x Temperature of the ambient fluid 

Tw Temperature at the surface 

{u,v) Dimensionless velocity components along the 
(x, y)axes 

V Wall suction velocity 7^-7^ 

(x, y) Axis in the direction along and normal to the 
surface respectively 

Greek symbols 

a Equal to % R d 

P Coefficient of thermal expansion 

A Equal to ^vv - ' 

AT Equal to 

'• Similarity variable 

6 Dimensionless temperature function 

{) w Surface temperature parameter 

p Viscosity of the fluid 

V Kinematic viscosity 

<5 Similarity variable 

p Density of the fluid 

a Stephman-Boltzman constant 

as Scattering co-efficient 

PS Absolute Viscosity at the film temperature 

T Coefficient of skin friction 

7w Shearing stress 

i// Non-dimensional stream function 

Subscripts 

w wall conditions 

Ambient temperature 

Introduction 

From technical point of view, the electrically 
conducting fluid flow in presence of magnetic field is 
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significant. For the last couple of decade, the 
electrically conducting fluid flow problems have 
received much attention from many researchers. The 
electrically conducting fluid over an surface that flows 
and transfers heat needs to study for the sake of 
resolving many mechanical, hydro dynamical 
engineering and astrophysical problems including 
solar structure, especially in the outer layers, the solar 
wind bathing and interstellar magnetic fields. 

The main application of magneto hydrodynamic fluid 
flow is in the astrophysical and geophysical field. 
However, as material walls are destroyed during 
creation and containment of hot plasmas during 
generating fusion power by electromagnetic forces, to 
resolve the problem, MHD fluid flow is applied. MHD 
free convection flow with variable viscosity from 
porous vertical plate has attracted attention not only 
for its fundamental aspects but also for its significance 
in the contexts of space technology involving high 
temperature. 

The interaction of the magnetic field and the moving 
electric charge carried by the flowing fluid induces a 
force, which tends to oppose the fluid motion. The 
magnetic force is proportional to the magnitude of the 
longitudinal velocity of fluid. Near the leading edge of 
the surface, the velocity is very small that the magnetic 
force acts in the opposite direction and is also 
minuscule. Consequently, the influence of the 
magnetic field on the boundary layer is exerted only 
through induced forces within the boundary layer 
itself. 

Merkin(1972) concluded free convection with blowing 
and suction. Lin & Yu(1988) studied free convection on 
a horizontal plate with blowing and suction. Hossain 
et al. [Hossain, Khanafer and Vafai (2001)] studied the 
effect of radiation on free convection flow with 
variable viscosity from a porous vertical plate. Hossain 
et al. [Hossain, Munir and Rees (2000)]studied flow of 
viscous incompressible fluid with temperature 
dependent viscosity and thermal conductivity past a 
permeable wedge with variable heat flux. Hossain 
&Takhar (1996) studied radiation effect on mixed 
convection along a vertical plate with uniform surface 
temperature. Molla et al. [Molla, Hossain and Yao 
(2004)]examined natural convection flow along a 
vertical wavy surface atuniform surface temperature in 
the presence of heat generation/ absorption. Akhter 
(2007) studied the effect of radiations on free 
convection flow on sphere with isothermal surface and 
uniform heat flux. Ali (2007)studied the effect of 
radiation on free convection flow on sphere with heat 


generation. Hossain et al. [Hossain, Alim and Rees 
(1999)] investigated the effect of radiation on free 
convection flow from a porous vertical plate. They 
[Hossain, Alim and Rees (1999)] performed a full 
numerical solution and found an increment in 
radiation parameter Rrfwhich makes the boundary 
layer thin and found an increment in surface 
temperature parameter which causes the boundary 
layer to be thick. Molla et al. [Molla, Hossain and 
Taher (2006)] studied the magnetohydrodynamic 
natural convection flow on a sphere with uniform heat 
flux in presence of heat generation. Gary et al. (1982) 
and Mehta and Sood (1992) concluded that when this 
effect is included, the flow characteristics substantially 
change compared to the constant case of viscosity. 
Recently, Kafoussius& Williams (1995) and 
Kafoussiaset al. [Kafoussias, Rees and Daskalakis 
(1998)] have investigated the effect of the temperature- 
dependent viscosity on the mixed convection flow past 
a vertical flat plate in the region near the leading edge 
using the local non-similarity method. In these studies, 
they concluded that when the viscosity of a fluid is 
sensitive to temperature variations, the effect of 
temperature-dependent viscosity has to be taken into 
consideration, otherwise considerable errors may 
occur in the characteristics of the heat transfer process. 
Hossain et al. [Hossain, Kabir and Rees (2002)] have 
investigated the natural convection of fluid with 
variable viscosity from a heated vertical wavy surface. 
Hossain& Munir(2000) investigated the mixed 
convection flow from a vertical flat plate for a 
temperature dependent viscosity. In the above studies 
[Hossain, Kabir and Rees (2002)] [Hossain & Munir 
(2000)] the viscosity of the fluid has been considered to 
be inversely proportional to a linear function of 
temperature. 

According to the state of the art, variable viscosity and 
MHD effects on laminar boundary layer flow of the 
fluids along porous plate has not been studied yet. The 
variable physical property of viscosity may change 
significantly with temperature. In this study the 
presence of variable viscosity on MHD free convection 
boundary layer flow from a porous vertical plate of a 
steady two dimensional viscous incompressible fluid 
has been considered. Here we have investigated the 
variable viscosity on MHD natural convection flow 
from a porous vertical plate numerically. The results 
are obtained for different values of relevant physical 
parameters and shown in graphs as well as in tables. 

The governing partial differential equations are 
reduced to locally non-similar partial differential 
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forms by adopting some appropriate transformations. 
The transformed boundary layer equations are solved 
numerically using implicit finite difference scheme 
together with the Keller Box technique[Keller (1978)]. 
Here, Our attention is focused on the evolution of the 
surface shear stress in terms of local skin friction and 
the rate of heat transfer in terms of local Nusselt 
number, velocity profiles as well as temperature 
profiles for selected values of parameters consisting of 
MHD M, variable viscosity y, Prandtl number Pr. 


Problem Description and Mathematical Model 

To investigate MHD free convection flow from a 
porous plate with variable viscosity, the fluid is 
assumed to be a grey, emitting and absorbing but non 
scattering medium and the surface temperature of the 
porous vertical plate, Tw, is constant, where T w > T m . 
The considered physical configuration is shown in Fig. 
1(A): 


x 



FIG. 1(A) THE COORDINATE SYSTEM AND THE PHYSICAL 
MODEL 


The conservation equations for the flow characterized 
with steady, laminar and two dimensional boundary 
layer, under the usual Boussinesq approximation, the 
continuity, momentum and energy equations can be 
written as: 


^ + ^=0 
dx dy 


(i) 


p (“ £ + » |) = t. O ' |) + psHt - r -> - ^° 2 “ < 2 > 

( dT dT\ , d z T 

P c p{ u r x + V ^) = k ^ & 

With the boundary conditions 

x = 0, y > 0, u = 0, T = T 00 - 

y=0, x>0, u=0, v=-V, T=T W 

y — » oo, x > 0, u = 0, T = T 
00 


where/? is the density, /3o is the strength of magnetic 
field, cm is the electrical conduction, k is the thermal 
conductivity, /? is the coefficient of thermal expansion, 
v is the reference kinematic viscosity v = /u! p, // is the 
viscosity of the fluid, C v is the specific heat at constant 
pressure. The absolute viscosity p is assumed to be 
vary with temperature according to a general 
functional form p = p : s(T), where pf is the absolute 
viscosity at the film temperature T/ and s(Tf) = 1. This 
form is chosen to allow definition of the stream 
function based on the absolute viscosity at the film 
temperature. For liquids, all transport properties vary 
with temperature. However, for many liquids like 
petroleum oils, glycerin, glycol, silicon fluids and some 
molten salt, the rate of variation of absolute viscosity 
with temperature is much more than that of the other 
properties. Under the above conditions, an analysis 
incorporating the above assumptions and describing 
the momentum and thermal transport within the flow 
field is more accurate than the usual assumption of 
constant properties evaluated at some reference 
temperature. It should be mentioned here that there 
are some fluids for which properties other than p vary 
strongly with temperature. In particular, water and 
methyl alcohol exhibit strong variation of both p and (i. 
The analysis presented here is not applicable to these 
liquids since only the variation of the absolute 
viscosity is taken into consideration as a function of 
temperature. However, for the case of an isothermal 
surface (in an un-stratified ambient fluid), the 
variation of the absolute viscosity with temperature 
takes the form p = luS( 0), where 6 is the dimensionless 
temperature in the boundary layer defined in equation 
(4), such that S(l/2) = 1. A wide variety of functional 
forms of S(9) satisfying this requirement was 
investigated in the literature such as algebraic 
expressions, power series, exponential forms, etc. 
Carey & Mollendorf (1978), the simplest form of the 
absolute viscosity is used in this investigation as 
follows: 


P~Pf 



(5a) 


This simple form amounts to a linear variation of the 
absolute viscosity with temperature, with the slope dp/ 
dT, evaluated at film temperature. The assumed linear 
variation of viscosity with temperature gives rise to a 
new parameter y defined by 



Now introduce the following non-dimensional 
variables: 
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hf e 2 


4x 

v-gfiKT 


' ,¥=V Vg/WT£ 3 j/+| 


f) ~^ w 


Where, 0 is the non-dimensional temperature function, 
Ow is the surface temperature parameter. 

Substituting (6) into Equations (1), (2) and (3) leads to 
the following non-dimensional equations 


l+/(0-q] f'" + 0-2f' 2 +3ff" + Zf" 
\ ^ J 


= VVWp v rr 

-1 (8) 
PrdT][dT}\ % d£ 84 ) I s ) 

Where Pr = pC P /k is the Prandtl number and 
M=j3o 2 cm/ vp is the hydromagnetic parameter. The 
boundary conditions (4) become 

f = 0, f = 0, 0 = 1 at t] = 0 
f = 0, 0 = 0 as r/ — > oo 

The solution of equations (6), (8) enable us to calculate 
the nondimensional velocity components u, v from 
the following expressions 

v -1 df (10) 

v = y=4 Qf + S-rjf’ + S-L) 

In practical applications, the physical quantities of 
principle interest are the shearing stress r«- and the rate 
of heat transfer in terms of the skin-friction coefficients 
C f x and Nusselt number Nux respectively, which can be 
written as 

NUx ~va 0 (ii) 

where r } and q =-*M 

W H= 0 W,=0 ( 12 ) 

q< is the conduction heat flux. 

Using the Equations (6) and the boundary condition (9) 
into (11) and (12), we get 


c fx =4+f)/'U°) 


Nu x =4 0'{xfi) 

The values of the velocity and temperature 
distribution are calculated respectively from the 
following relations: 

«=£ 2 fit'll 0=0(x,y) ( 14 ) 


Numerical Procedure 

Using implicit finite difference method with Keller-Box 
Scheme, the solution of the local non similar partial 
differential equation (7) to (8) subjected to the 
boundary condition (9) are obtained, which has been 
described in details by Cebeci and Bradshaw(1991). 

The solution methodology of equations (7) and (8) 
with the boundary condition given in equation (9) for 
the entire £ values based on Keller-box scheme is 
proposed here. The scheme specifically incorporated a 
nodal distribution favoring the vicinity of the plate, 
enabling accuracy to be maintained in this region of 
steep gradient. In detail, equations (7) and (8) are 
solved as a set of five simultaneous equations. 


!+rf + 3 #" - 2 (f) 2 + 0~ If" ~ Mf’4 2 


= £ r §rjZ. r 

h { 84 84 J 


ilM + ^ e 'A r wfc) (16) 

To apply the aforementioned method, we first convert 
Equations (15)-(16) into the following system of first 
order equations with dependent variables m (4t , 7) 
v(#,7 ), Pi£,V) and g{£,n) as 

/= u, u'= v, g=6, and 9 = p (17) 


CO 2 2 

1+ /> 5 g-y v +p l fv-p 2 u +g-4v-p 4 u4 


du df 

l — 1 




where 


pi =3, p 2 = 2, pt = M and ps = y (20) 

The corresponding boundary conditions are 

f(4, 0) = 0, u(4, 0) = Oandg(£, 0) = 0 

u(4, oo) = 0,g(£,oo) = 0 ( 2 i) 

We now consider the net rectangle on the ( c, >]) plane 
and denote the net point by 

, ? 0 =0 ’ r,j=r,j_ x + h j, j = 1 , 2, ...7 

4° = 0 , 4 n = 4 n ~ l + k , n= 1,2, N 

Here ‘n’ and ‘j’ are just sequence of numbers on the 
(<fj, rj) plane, kn and hj are the variable mesh widths. 

We approximate the quantities (f, u, v, p) at the points 
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(£" > Vj ) of the net by (•// > 11 j » v j > Pj ) which we call net 
function. 


1 


s n - yz =\(s n +e n - 1 ) 

*/~ 1/2 = + -V" 1 ) 

«y - 1/2 " W + *3-i } 


( 22 ) 



FIG. 1(B) NET RECTANGLE FOR DIFFERENCE 
APPROXIMATIONS FOR THE BOX SCHEME 

Now we write the difference equations that are to 
approximate Equations (17)-(19) by considering one 

mesh rectangle for the mid-point (£ > n j- \C) to obtain 


f n - f n 

S i J J- ] - u n 

~ u j-\l2 


h j 


u”-u” , 
■/ i " 1 

A . 

J 


g n -g' 7 , 

~ p j- 1/2 


(23) 

(24) 

(25) 


7 


Similarly, Equations (18)— (19) are approximate by 


centering about the midpoint 7 ^ » "J / 2 


7 2 , « /_ 1 / 

Centering the Equations (22) about the point 


(£ W /2,n) without specifying 77 to obtain the 
algebraic equations. The difference approximation to 
Equations (18)-(19) becomes 

'' _1 {l+P 5 (g-0.5)v}"-{ 1 +p 5 (g- 0 .5)v}"_ i 

+{ 1 + P 5 (g-0.5)v}" 1 -{l+P 5 (g-0.5)v}"_[ 

+ «P \ )"_ 1 / + % 1 / - «/» 2 1 / + «» K“ 2 )"_ 1 / 

7 72 7 72 7 /2 7 /2 


1 / 


7 

2 


1 n — 1 

72 2 2 n n «— 1 

hi *“" ! 7-xs-x 


72—1 


72—1 


7 /2 7 72 7 72 


where 

72—1 -1 

I ! =*■ 

'-2 


{l+p 5 ( g-0.5 )}" ' - -{!+P 5 ( . *-0.5 )}"_j 


J ~2 ] ~2 


n - 1 




7 2 7 2 7 2 


+*y \fj )j_ x 

J 2 

and 


-1 


,72— 1 1 


2 72-1 


7? 


_ n-1 n - 1 

-L j + a #) 

7 2 7 2 


7 2 7-- 

,3 ,n 


=> ^ [*y 1 (P” - ) + h j 1 { + A g) 3 l" 

-{p 3 p(l + Ag) 3 }"_!}] + ^" wp” ] / + (Px ) n w(fP )” l7 

7 7 -/2 7 -/7 1 7-/7 7"7 

!— 1 n — 1 /? — 1 

1 / + an [—(wg) | / + ifp) | / ] 

72 7 72 7 72 


n — i 

= -M , 

J-/2 


+a 


Mug)” x/Hfpf. v -u n 1/^“ 1 1/+“"' 1 1/ 

7_ /7 7_ /7 7 /2 7_ /7 7 /7 

tj -1 n 


77—1 , 77-1 n 1 

•+«. i/g. 1 /} 


7-72 


77 77-1 77-1 77 

+ p j-yfj-y ~ p j-yfj.y} 

7 /2 7 /2 7 /2 7 72 


1—1 n n n n 2 n n 

PF [*y (Py - Py_i )1 + # , 1/ P , 1/ + (P 4 f ) 1/ * 1/ 

7 / 2 7 / 2 7 / 2 7 / 2 


1 -1 
pF^ 

+{(Pl)"_ I/^H/P)"!/ 

7 72 7 72 

^ rf/ \72 / \72— 1 72 72—1 . 72- 1 72 i 

[{(»g) , / -(«g) i/ -« . i/ g . i/ +« . i/ g . i/ 1 

7 72 7 72 7 72 7 72 7 72 7 72 

7 72 7 72 7 72 7 72 

«-l 

= T y 

7 72 

where 

«— 1 

M W 

j A 


-Pj- \ 


1 —1 /I— 1 «-l 

IT ; '/ (/ ’/ /’z 

-N"i/ p'^z+cp/r 1 !/ (/p)”i/ ] 

7_ /2 7_ /2 7 -/^ 7_ /2 

77—1 77-1 72—1 

1 1/ = - M 1/ + K/ 1/ -<“«) J/ 1 

7 72 7 72 7 72 7 72 

The corresponding boundary conditions (21) become 

/« = 0 , 77» = 0 , g» = 1 

zz} = o, g} = 0 
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which just express the requirement for the boundary 
conditions to remain during the iteration process. Now 
the momentum and energy equations will be 
converted into system of linear Equations and together 
with the boundary conditions can be written in matrix 
or vector form, where the coefficient matrix has a block 
tri-diagonal structure. The whole procedure, namely 
reduction to first order followed by central difference 
approximations, Newton's quasi-linearization method 
and the block Thomas algorithm, is well known as the 
Keller-box method. 

Results and Discussion 

To investigate the MHD natural convection flow from 
a porous vertical plate with variable viscosity, 
numerical values of local rate of heat transfer was 
calculated in terms of Nusselt number Ntix for the 
surface of the porous vertical plate from lower 
stagnation point to upper stagnation point. Simulation 
was performed for different values of the 
aforementioned parameters and those are shown in 
tabular form in Table 1 and graphically in Figure 5-7. 
The effect for different values viscosity yor\ local skin 
friction coefficient Ch and the local Nusselt number 
Nux, as well as velocity and temperature profiles are 
displayed in Figure 2 to 7. 



FIG. 2(a) VELOCITY AND (b) TEMPERATURE PROFILES FOR 
DIFFERENT VALUES OF VISCOSITY PARAMETER T WITH 
OTHERS FIXED PARAMETERS. 

Keeping constant the Prandtl number Pr = 1.0, 
magnetohydrodynamic parameter M = 1.0 surface 
temperature parameter 0w = 1.1 and figures 2 (a)- 2 (b) 
display results for the velocity and temperature 
profiles for different values of viscosity parameter y= - 
1.0, 0.0, 0.5, 1.0, 1.9. Figures 2(a)-2(b) shows that as the 
viscosity parameter y increases, the velocity profiles 
decreases and the temperature profiles increase. The 
fact is that due to the viscosity it takes away the warm 
fluid and thereby decreasing the maximum velocity 
with a decreasing intensity of the natural convection 
rate as depicted be figure 2 (a). The velocity is zero at 
the boundary wall then the velocity starts increasing to 
the peak value while r| increases from 1 to 1.8 and then 
starts decreasing after 77 = 1.8 and so on. Finally the 
velocity approaches to zero (the asymptotic value). 




FIG. 3(a) VELOCITY AND (b) TEMPERATURE PROFILES FOR 
DIFFERENT VALUES OF PRANDTL NUMBER PR WITH OTHERS 
FIXED PARAMETERS. 

However, in figures 3(a)-3(b) it has been shown that 
keeping constant 6L = 1.0, M = 1.0 and y= 0.5 while the 
Prandtl number Pr (= 0.7, 1.0, 2.0, 3.0 and 4.0)increases 
both the velocity and temperature profiles decrease. 
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This is due to the fact that as Pr increases, viscosity 
increases so the velocity decreases. 

Figure 4(a) shows the velocity profiles for varying M 
from 0 to 30 with Prandtl number Pr = 1.0, viscosity 
parameter y= 1.0 and surface temperature parameter 
(L = 1.1. It is shown in figure 4(a) that as the MHD 
parameter increases the velocity profiles decrease 
because electrically conducting fluid affects the flow so 
the velocity decreases. It is also observed from figure 
4(a) that the changes of velocity profiles in the r| 
direction reveal the typical velocity profile for natural 
convection boundary layer flow. With an increment 
ofr), the velocity profile shows an asymptotic 
behaviour. The maximum values of velocity are 
recorded to be 0.11934 and 0.14053 at p = 0.88811 and 
0.16893, 0.18645, 0.20630 at ?/= 0.94233 for M =.30.0, 
20.0,10.0, 5.0 and 0.0 respectively. The velocity is 
0.20630 at r/= 0.94233 for M = 0.0. Here, it is observed 
that at rp 0.88811, the velocity decreases by 73.40% as 
the MHD parameter M changes from 0 to 30.0. 




FIG. 4(a) VELOCITY AND (b) TEMPERATURE PROFILES FOR 
DIFFERENT VALUES OF MHD PARAMETER M WITH OTHERS 
FIXED PARAMETERS. 

Figure 4(b) shows results for the temperature profiles, 
for different values of MHD parameter M while 


Prandtl number Pr = 1.0, viscosity parameter y = 1.0 
and surface temperature parameter 61. = 1.1. From 
figure 4(b) it infers that as the MHD parameter M 
increases, the temperature profiles increase. It is 
observed that the temperature profile decreases 
gradually along ? 7 direction from 1.0 to 0 with an 
asymptotic behaviour. For M = 30.0, 20.0, 10.0, 5.0, 0.0 
the temperature profile raise to 0.52579, 0.50164, 
0.47181, 0.45488 and 0.43689 at r]= 0.88811 then it starts 
decreasing. And for other values, thoseare gradually 
increasing. 




FIG. 5 (a) SKIN FRICTION AND (b) RATE OF HEAT TRANSFER 
FOR DIFFERENT VALUES OF VISCOSITY PARAMETER T WITH 
OTHERS FIXED PARAMETERS. 

Figure 5 (a) shows that skin friction coefficient C/x 
increases with viscosity parameter considering 
Prandtl number Pr = 1.0, surface temperature 

parameter Qw = 1.1 and MHD parameter M = 1.0. It is 
observed from figure 5(a) that the skin friction 
increases gradually from zero at lower stagnation 
point along the c direction. Figure 5(b) reveals that the 
rate of heat transfer decreases along the £ direction for 
y= -1.0, 0.0, 0.5, 1.0 and 1.9. A hot fluid layer is created 
adjacent to the interface of the wall due to the viscosity 
mechanism and ultimately the resultant temperature 
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of the fluid exceeds the surface temperature. 
Accordingly, the heat transfer rate from the surface 
decreases as shown in Fig. 5(b). 

The variation of the local skin friction coefficient Cfx 
and local rate of heat transfer Nux for different values 

of Prandtl number Pr while 0 W = 1.0, M = 1.0 andy= 1.0 
are shown in Figures 6(a)-6(b). It can be observed from 
these figures that as the Prandtl number Pr increases, 
the skin friction coefficient decreases and the rate of 
heat transfer increases. 



Z, 



FIG. 6(a) SKIN FRICTION AND (b) RATE OF HEAT TRANSFER 
FOR DIFFERENT VALUES OF PRANDTL NUMBER PR WITH 
OTHERS FIXED PARAMETERS . 

Figures 7(a)-7(b) show that skin friction coefficient Cfx 
and heat transfer coefficient Nux decrease for 
increasing values of MHD parameter M while 
viscosity parameter y= 1.0, Prandtl number Pr = 1.0, 
and surface temperature parameter dw = 1.1. The 
values of skin friction coefficient Cfx and Nusselt 
number Nu x are recorded to be 0.25646, 0.30238, 
0.39508, 0.50394. 0.87147 and 1.00382, 1.00463, 1.00414, 
1.00463, 1.04285 for M = 30.0, 20.0, 10.0.5.0, 0.0 
respectively which occur at the same point £ = 1.5. 
Here, it observed that at £ = 1.5, the skin friction 
decreases by 70.57% and Nusselt number Nux 


decreases by 3.74% as the MHD parameter M changes 
from 0.0 to 30.0. It is observed from figure 7(a) that the 
skin friction increases gradually from zero at lower 
stagnation point along the /(direction. From Figure 7(b) 
it is revealed that the rate of heat transfer decreases 
along the £, direction. 



FIG. 7(a) SKIN FRICTION AND (b) RATE OF HEAT TRANSFER 
FOR DIFFERENT VALUES OFMHD PARAMETER Q WITH 
OTHERS FIXED PARAMETERS 

Numerical values of rate of heat transfer Nux and skin 
friction coefficient Cfx are calculated from Equations 
(13) from the surface of the vertical porous plate. 
Numerical values of C/vand Nux are shown in Table 1. 


TABLE 1 SKIN FRICTION COEFFICIENT AND RATE OF HEAT TRANSFER 
AGAINST \ FOR DIFFERENT VALUES OF MHD PARAMETER M WITH OTHER 
CONTROLLING PARAMETERS Pr = 1.0, Gw =1.1. AND y = 1.0. 


5 

M = 

= 30.0 

M = 

= 20.0 

Cfx 

Niix 

Cfx 

Nux 

0.01 

0.00720 

55.24526 

0.00720 

55.25031 

0.05 

0.03579 

11.37380 

0.03586 

11.38889 

0.10 

0.07033 

5.86455 

0.07081 

5.89219 

0.50 

0.23492 

1.34271 

0.26201 

1.41711 

1.00 

0.25731 

1.01311 

0.30335 

1.01638 

1.50 

0.25646 

1.00382 

0.30238 

1.00463 

5 

M = 

= 10.0 

M 

= 0.0 

Cfx 

Nux 

Cfx 

Nux 

0.01 

0.00720 

55.2553 7 

0.00720 

55.26042 

0.05 

0.03593 

11.40399 

0.03600 

11.41913 

0.10 

0.07130 

5.92000 

0.07180 

5.94800 
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0.50 

0.29917 

1.51681 

0.35101 

1.64306 

1.00 

0.39248 

1.03533 

0.65607 

1.16318 

1.50 

0.39508 

1.00414 

0.87147 

1.04285 


In the above table, the values of skin friction coefficient 
Cfx and Nusselt number Nux are recorded to be 0.25646, 
0.30238, 0.39508, 0.87147 and 1.00382, 1.00463, 1.00414, 
and 1.04285 for M = 30.0, 20.0, 10.0, 5.0, 0.0 and 
respectively which occur at the same point £ = 1.5. 
Here, it is observed that at £ = 1.5, the skin friction 
decreases by 70.57% and Nusselt number Nux 
decreases by 3.74% as the MHD parameter M changes 
from 0.0 to 30. 

Comparison of the Results 

In order to verify the accuracy of the present work, the 
values of Nusselt number and skin friction for Rj= 0.05, 
Pr = 1.0, y = 0, M = 0 and various surface temperature 

parameter 0 W = 1.1, 0 W = 2.5 at different position of 
are compared with Hossain et al. (Hossain, Alim, and 
Rees) as presented in Table 2. The results are found to 
be in excellent agreement. 


TABLE 2 COMPARISON OF THE PRESENT TABLE WITH HOSSAIN ET 
AL(Hossain, Alim, and Rees) 



0w= 1.1 

# 

Hossain 

Present [ 


Cfx 

Nux 

Cfx 

Nux 

0.1 

0.0655 

6.4627 

0.06535 

6.48306 

0.2 

0.1316 

3.4928 

0.13138 

3.50282 

0.4 

0.2647 

2.0229 

0.26408 

2.03018 

0.6 

0.3963 

1.5439 

0.39519 

1.55522 

0.8 

0.5235 

1.3247 

0.52166 

1.32959 

1.0 

0.6429 

1.1995 

0.64024 

1.20347 

1.5 

0.8874 

1.0574 

0.88192 

1.06109 


(hv = 2.5 | 

£ 

Hossain 

Present j 


Cfx 

Nux 

Cfic 

Ntlx 

0.1 

0.0709 

8.0844 

0.07078 

8.10360 

0.2 

0.1433 

4.2858 

0.14313 

4.29682 

0.4 

0.2917 

2.4003 

0.29120 

2.40669 

0.6 

0.4423 

1.7863 

0.44145 

1.78912 

0.8 

0.5922 

1.4860 

0.59080 

1.48991 

1.0 

0.7379 

1.1098 

0.73590 

1.31822 

1.5 

1.0613 

1.1098 

1.05693 

1.11262 


Conclusion 

For different values of relevant physical parameters 
including the viscosity parameter y MHD natural 
convection with variable viscosity flow from a porous 
vertical plate has been investigated and the following 
conclusions may be drawn: 

• Significant effects of MHD parameter M and 
viscosity parameter y on velocity and 
temperature profiles as well as on skin friction 


coefficientC/vand the rate of heat transfer Nux 
have been found in this investigation but the 
effect of MHD parameter M and viscosity 
parameter y on rate of heat transfer is more 
significant. An increase in the values of 
viscosity parameter y leads to the decrement of 
velocity and the increment of temperature 
profiles and the local skin friction coefficient Cfx 
increase and the local rate of heat transfer Nu x 
decreasing at different position of b' f o r Pr =1.0. 

• For increasing values of Prandtl number Pr 
leads to decrease of the velocity profile, the 
temperature profile and the local skin friction 
coefficient Cfx but increment of the local rate of 
heat transfer Nux. 

• An increase in the values of M leads to increase 
of the temperature profiles and the velocity 
profiles, the local skin friction coefficient Cfx 
and the decrease of local rate of heat transfer 
Nux. 
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